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We describe mechanical metamaterials created by folding flat sheets in the tradition of origami, the art of 
paper folding, and study them in terms of their basic geometric and stiffness properties, as well as load 
bearing capability. A periodic Miura-ori pattern and a non-periodic Ron Resch pattern were studied. 
Unexceptional coexistence of positive and negative Poisson's ratio was reported for Miura-ori pattern, 
which are consistent with the interesting shear behavior and infinity bulk modulus of the same pattern. 
Unusually strong load bearing capability of the Ron Resch pattern was found and attributed to the unique 
way of folding. This work paves the way to the study of intriguing properties of origami structures as 
mechanical metamaterials. 

Mechanical metamaterials, the man-made materials with mechanical properties mainly defined by their 
structures instead of the properties of each component, recently have attracted great attention 1 \ 
Origami, creating three-dimensional (3D) structures from two-dimensional (2D) sheets through a 
process of folding along creases, provides an interesting source for designing mechanical metamaterials and 
has been transformed by mathematicians, scientists, and engineers to utilize the folded objects' deformability and 
compactness in applications ranging from space exploration (e.g., a foldable telescope lens 5 ), to automotive safety 
(e.g., airbags 6 ), biomedical devices (e.g., heart stent 7 ), and extremely foldable and stretchable electronics 8,5 . 
Notable progress has been made in the area of origami theory including tree theory 10 and its corresponding 
computer program 11 , folding along creases 1214 , and geometric mechanics of a periodic origami pattern 15 . Among 
classes of origami patterns, a particular one, namely rigid origami, in which the faces between the creases remain 
rigid during folding/unfolding and only the creases deform, is different from most origami patterns that require 
faces bending or partial crumpling to make many-step folds. Idealized rigid origami possesses one of the most 
obvious advantages of origami in terms of deformation, i.e., the deformation is completely realized by the folding/ 
unfolding at the creases and does not involve any deformation at the rigid faces 4 . The geometric characteristics, 
such as the necessary condition around a single vertex for rigid origami 1617 have been studied, and a computer 
simulator for rigid origami 18 exists. There have also been made limited efforts to study the structural character- 
istics of one particular rigid origami, namely Miura-ori 19 , as a mechanical metamaterial, with the main focus on 
the negative Poisson's ratio 15-20 , though these properties can be more rigorously examined. It is noticed that the 
existing studies are mainly focused on periodic origami patterns (e.g., Miura-ori); however, non-periodic origami 
patterns as mechanical metamaterials have not gained attention yet, partially due to the difficulties in theoretical 
analysis and modeling. To span a much wider spectrum of using rigid origami as mechanical metamaterials, we 
report a systematic study of two rigid origami folding patterns, not only the periodic Miura-ori but also a non- 
periodic Ron Resch folding 21 using combined analytical and numerical approaches. Specifically, we rigorously 
address the commonly mistaken in-plane Poisson's ratio of Miura-ori, which was believed to be always negative 
but we show here that it can be positive as well, and its physical interpretation. The ubiquitous non-local 
interactions between vertices of rigid origami patterns are captured through a non-local finite element approach 
and the compressive buckling resistance of a Ron Resch tube is studied, which inspires a theoretical and 
experimental study of the load bearing capability of the Ron Resch pattern. The result witnesses the superb load 
bearing capability of this Ron Resch pattern. Based on the approaches in this paper, mechanical properties of 
different rigid origami patterns, both periodic and non-periodic ones, can be readily studied. 

Results 

Unit cell and the whole pattern of a Miura-ori. Figure la illustrates a Miura-ori (n lt n 2 ) in its folded state, 
containing n x ( = 11) vertices in x t direction and n 2 ( = 11) vertices in x 2 direction, with % 3 as the out-of-plane 
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(a) 



"Mountain" crease 
"Valley" crease 



(b) 





(c) 



Central vertex 
(d) 



Figure 1 | Illustrations of Miura-ori. (a) A Miura-ori (rii, n 2 ) in its folded state with n x vertices in X\ direction, n 2 vertices in x 2 direction. x 3 is the 
out-of-plane direction. Specifically for this illustration, tii = 11, » 2 = H> ft = 45° and a/b= 1/V2. (b) A Miura-ori in its planar state, corresponding to 
(a). The solid lines represent "mountain" creases that remain on the top after folding. The dashed lines represent "valley" creases that remain on the 
bottom after folding, (c) A unit cell of a Miura-ori. a 1 ,a 2 are two dihedral angles. In each parallelogram, the length of the short side is a and that of the long 
side is b, with the acute angle of /?. The projected angle between the two ridges is <j>. The size of the unit cell is /, w, and ft, in X\, x 2 , and x 3 directions, 
respectively, (d) A non-local element for Miura-ori that focuses on the central vertex. 



direction. Its corresponding planar state is shown in Fig. lb. The 
geometry of a Miura-ori is defined by many identical rigid 
parallelogram faces (with four gray ones highlighted in Fig. lb) 
linked by edges that can be folded into "mountain" and "valley" 
creases. The Miura-ori is a periodic structure and its unit cell is 
shown in Fig. lc, where the four parallelograms are identical with 
the short sides of length a, the long sides of length b, and the acute 
angle [i £ [0°,90°] . Since the necessary condition for rigid origami 16,17 
states that there are n — 3 degrees of freedom, where n is the number 
of edges at one vertex, Miura-ori with n = 4 has only one degree of 
freedom. Therefore, if the shape of a parallelogram face is prescribed, 
i.e. p, a and b are given, one parameter <j>e [0°,2/J], the projection 
angle between two ridges, can be used to characterize the folding of 
the unit cell of Miura-ori, with (j) = 2/i for the planar state and (f> = 0° 
for the completely collapsed state. The size of the unit cell is I = 



and equal 180° for the planar state and 0° for the completely 
collapsed state. When the whole structure of a Miura-ori is put in 
an imaginary box with the dashed lines as the boundaries (Fig. la), 
the dimension of the whole Miura-ori is then given by 



L=(«i-l)bsin0/2), 
cos p 



W={n 2 -l)i 



s^/2) 



+ bcos{<f>/2), 



(2) 



ai/sin 2 /?— sin 2 (</;/2) 
cos(^>/2) 



cos P 



i\J sin 2 P~ 



2bsm(<p/2), w = 2a i-rj-z, and h= t— — , in x lt 

cos(0/2) cos(<p/2) 

x 2 , and x 3 directions, respectively. It is noted that the length of the 

"tail" bcos(<[>/2) is not considered in the unit cell 15 . The periodicity of 

this pattern only requires two dihedral angles a 1 S [0°,180°] and 

a 2 G[0°,180°] to characterize the geometry (Fig. lc), which are 

given by 



V-l = cos 



1-2- 



(1) 



a 2 = cos 1 [l-2 cot 1 p tan 2 (<j>/2) 



and thus the imaginary volume occupied by this Miura-ori is given by 
V = LxHxW. (3) 



Apparently even the Miura-ori is periodic, its size in x 2 direction 
(i.e. W) is not proportional to its counterpart for the unit cell, w, due 
to the existence of the "tail" with length fccos(^i/2). Consequently, it is 
not accurate to use the unit cell to study the size change of a whole 
Miura-ori (e.g., Poisson's ratio), particularly for smaller patterns 
(e.g., used in 15,20 ). 

In-plane Poisson's ratio of Miura-ori. In-plane Poisson's ratio of 
Miura-ori is believed to be negative from intuitive observations and 
as testified by some theoretical studies using the unit cell (Fig. lc) 15,20 . 
An accurate mean to define the Poisson's ratio is to use the size of a 
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whole Miura-ori, instead of using the unit cell. Specifically, the in- 
fill 



plane Poisson's ratio v 12 is defined as v 12 =- 



<' : 22 



where 



822 ->o 



dL dw 

£11 = — and s 2 2 = are the infinitesimal strains in X\ - and x 2 - 

L W 

directions, respectively. Using equation (2), the in-plane Poisson's 
ratio v 12 is obtained as 



cot 2 (^/2) 



(n 2 — cos/?+ cos 2 (</>/2) 
(n 2 — l)f/cos/?— cos 2 ((f>/2) 



(4) 



where (7 = alb. Another in-plane Poisson's ratio v 2 \ is just the 
reciprocal of v 12 . Figure 2a shows the contour of v 12 as a function 
of angle <j) and a combination parameter (« 2 — l)!7cos/?. Clearly, v 12 
can be negative or positive, which is different from commonly 




observed negative in-plane Poisson's ratio. The boundary 
separating the negative and positive regimes of v 12 is defined by 
vanishing the denominator of v 12 , i.e., (« 2 — l^fjcos/? = cos 2 (</>/2). 
At this boundary, v 12 flips between positive and negative infinities; 
thus V12 e [ — co , + cc ]. For one scenario, where n 2 = 5 (small pattern), 
11 = 1/2, [3 = 78.5°, and thus (n 2 — l))7cosjS < 1, v 12 is positive for 
i^e [0,101.5°] and changes to negative for ^ e [101.5°,2/?]. For 
another scenario, « 2 = 13 (large pattern), 17 = \/\fl, ji = 45°, and 
thus (n 2 — l)f7cos/? > 1, v 12 stays negative, as reported by others 
using the unit cell 15,20 ; and the Miura-ori becomes an auxetic 
material. Similar analysis can be applied on the out-of-plane 
Poisson's ratios. Details can be found in the Supplementary 
Information, Section "Out-of-plane Poisson's Ratio". 

Figure 2b provides an intuitive explanation of the sign change in 
the in-plane Poisson's ratio v 12 . For the specific example with rii = n 2 
= 5, 1] = 1/2 and /3 = 78.5°, the size of this Miura pattern in X\ 
direction, L, decreases monotonically from the planar state to the 
collapsed state, which is pictorially shown in the three insets for <j> = 
157°( = 2/i), (j) = 140°, and <j) = 20° with U> L 2 > L 3 . In contrast to 
L, the respective size of this pattern in the x 2 direction, W, does not 
change monotonically with the angle tf>. From the planar state to the 
collapsed state, Wi > W 2 but W 2 < W 3 , which gives v 12 < 0 when L 
and W change in the same direction and v 12 > 0 when L and W 
change in the opposite direction. The non-monotonic change of Wis 
due to the "tail" term bcos((j>/2) in equation (2), which was missed in 
previous studies 15,20 . As shown in the Supplementary Information, 
Section "Change of Length W" for more details, the two terms in W 
(equation (2)) dominate at different stage of folding. 

In addition to the negative and positive in-plane Poisson's ratio of 
the Miura-ori, the ranges of Poisson's ratios, specifically, 
v 12 e [— °o,+=°], v 13 e[0,°°], v 23 e[— °c,=°] (see Supplementary 
Information, Section "Out-of-plane Poisson's Ratio") are also fas- 
cinating if the range of Poisson's ratio for common materials is 
considered as the reference, i.e., v6[- 1,0.5] . Now we interpret these 
fascinating phenomena in terms of shear and bulk modulus of 
Miura-ori. 



-,4 




Ll i. 



(b) 

Figure 2 | Poisson's ratios of Miura-ori. (a) Contour plot of in-plane 
Poisson's ratio v 12 as a function of (j) and the combined parameter 
(n 2 — 1)j7Cos/?. (b) Explanation of negative and positive in-plane 
Poisson's ratio v 12 . 



Miura-ori subjected to shear and hydrostatic deformation. To 

study the shear deformation that is non-uniform across the Miura- 
ori, we developed a numerical approach to characterize the geometric 
features of the Miura-ori, i.e., the non-local interactions between 
rigid faces. As shown in Fig. lb, the vertex marked by the solid 
blue dot not only interacts with its nearest-neighboring vertices 
(marked by the solid red dots) through the rigid faces, but also its 
second-neighboring vertices (marked by the solid green dots) 
through dihedral angles. Thus the interactions between vertices are 
non-local. This non-local nature is ubiquitous in rigid origami and 
can be more complicated for other patterns, which can be further 
illustrated by the Ron Resch pattern 21 , detailed in the Supplementary 
Information, Section "Non-local Interactions in the Ron Resch 
Pattern" (e.g., Supplementary Fig. S3). We developed a non-local 
finite element based model and Fig. Id shows the non-local 
element for Miura-ori. Details can be found in Methods. 

Supplementary Figure S4 shows a deformed state of a (n l = 13, n 2 
= 13) Miura-ori subjected to a finite shear force in the negative Xi 
direction. Here it is noticed that an initially periodic Miura-ori 
deforms non-uniformly under shear loading, which disables the def- 
inition of a shear modulus. It is seen that the Miura-ori responds in 
an opposite way to shear force. Specifically and clearly, the vertical 
lines tilt to the positive Xi direction, as the shear force is applied along 
the negative Xi direction. This opposite relationship is thus consistent 
with v 12 e [—oo ; +00]. 

The bulk modulus K of Miura-ori can be defined the same way as 
that in continuum mechanics to link the hydrostatic pressure p and 
the volumetric strain 0( = e u + £ 22 + £ 33 ), 
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(5) 



p = 0 



Using the principle of superposition (details given in the Supple- 
mentary Information, Section "Bulk Modulus of Miura-ori"), the 
bulk modulus K is given by 

1 l-v„-v 31 l-Vi2-v 32 I-V13-V23 

1 ~ 1 s , (°j 



K 



where En 



£11 



V12- 

-E22 



•V13- 
£33 



-B22 — 



dff2 



da 33 
de 



are the 



ife, 3 =0 



<feii=0 Mt 22 ^ = 0 "^33 

tangential moduli of the stress-strain curve. Using the work conjug- 

„ , , SW t ot dW tot 
ate relation, stresses are expressed as <Tn=— , ct 2 2 = 



SW tc 



8s n ' ~" <3£ 22 

(T33 = """ , where W tot = U tot /Vis the elastic energy density with 

c% 3 

l/ to( given by Supplementary Eq. S2 and V given by equation (3). As 
shown in the Supplementary Information, Section "Range of Tensile 
and Bulk Modulus" for details, the tensile (-En, £ 2 2> £33) an d bulk 
moduli (K) have a wide range of variation and some of them vary 
from 0 to infinity, such as K. 

Ron Resch pattern and its buckling resistance. Next we study a 
non-periodic rigid origami folding, namely a Ron Resch pattern, 
using the developed non-local finite element approach. The Ron 
Resch pattern and its non-local elements are given in 
Supplementary Fig. S3. To illustrate the non-periodicity, several 
Ron Resch patterns (specifically, a Ron Resch dome, a tube and a 
stingray) have been studied and the histograms of the three dihedral 
angles ft, ft, and ft are shown in Supplementary Fig. S5. It is 
obvious that the Ron Resch pattern is non-periodic and the 
importance of a universal numerical platform to study this type of 
rigid origami is thus apparent. We first study the buckling resistance 
of a Ron Resch tube (Fig. 3a). A Ron Resch tube in its folded state 
contains many equilateral triangles. As shown in the zoom- in details 
in Fig. 3a, the dihedral angles ft G [0°,90°] and ft G [90°,180°]. 
Because of the folded state, the centroids of these equilateral 
triangles form spikes pointing to the central axis of the Ron Resch 
tube as shown in the top view of Fig. 3a. The boundary condition for 
the axial compressive buckling is that one end of the tube is fixed and 
the other is subjected to a compressive force, which is the same as the 
Euler buckling. Figure 3b shows the compressive force normalized by 
k RR /b varies as the compressive strain increases, and the insets show 
some characteristic snapshots at the compressive strains of 13%, 30% 
and 45% from left to right, respectively. Here k m is the spring 
constant of the hinges for dihedral angles (detailed in the 
Supplementary Information, Section "Work Conjugate Relation - 
Stress and Moduli for Miura-ori"), and b is the size of the right 
triangles (Supplementary Fig. S3). It is interesting to find that 
buckling does not occur, which can be explained by the negative 
Poisson's ratio. Upon compression, the two dihedral angles ft and 
ft decrease, which lead to the further pushing the spikes towards the 
central axis of the tube. Thus the compression leads to a shorter tube 
with smaller radius due to the negative Poisson's ratio (the leftmost 
inset of Fig. 3b). Further compression leads to an even smaller tube 
radius (the middle inset of Fig. 3b). Eventually, the equilateral 
triangles form completely folded states, which is captured by ft = 
0°, ft = 120° and results in a much smaller tube radius (the 
rightmost inset of Fig. 3b). At the completely folded state, the tube 
cannot be further compressed because of the rigidity. Thus, axial 
compressive force does not lead to the buckling of a Ron Resch tube. 

Load bearing capability of a Ron Resch plate. This intriguing 
buckling resistance phenomenon motivates a further study of the 



load bearing capability of the Ron Resch pattern. The compressive 
load applied on top of a Ron Resch dome leads to a completely 
compact and flat state (namely, a Ron Resch plate), where the 
equilateral triangles collapse to three-fold structures with ft = 0°, 
ft = 120°, and ft = 90° (Fig. 3c). Figure 3d shows the striking load 
bearing capability of a Ron Resch dome folded from a single sheet of 
20-lb copy paper: a 32.4 lb load is carried by a Ron Resch plate with 
actual mass 4.54 g. This remarkable capability is mainly a result of 
the folded structure, not the material properties of the paper, which 
suggests that origami can produce exceptional mechanical 
metamaterials. Figure 3e shows snapshots of the bottom of Ron 
Resch plate when 3Tb load (left panel) and 32.4-lb load are applied 
(right panel). It is found that at the failure point (where the 32.4-lb 
load is applied), the tips of the three-fold structures are flattened and 
instability occurs. To compare the load bearing capability of a Ron 
Resch plate with three-fold supporting structures and commonly 
seen six-fold ones that are used in airplane wings, the buckling 
analysis is conducted to compute their critical compressive loads 
P cr by using finite element package ABAQUS (details are given in 
the Supplementary Information, Section "Buckling Analysis of a Ron 
Resch Plate and a Six-Fold Supporting Structure"). Figures 3f and 3g 
show the first buckling modes for the Ron Resch plate and the six- 
fold structure. By assigning the same geometric parameters 
(including thickness and height of the support) and material 
properties (including elastic modulus and Poisson's ratio), P cr of 
the Ron Resch plate is about 50% larger than that of the six-fold 
structure. Though the six-fold structure has higher symmetry to 
increase P cr in a linear fashion (i.e., P cr ~ order of symmetry), the 
decreasing height of the support for the Ron Resch plate from the 
center to the surroundings increases P cr in a quadratic fashion (i.e., 
P cr ~ 1/height 2 ), which endows a higher load bearing capability of the 
Ron Resch plate. This result suggests that generalized Ron Resch 
patterns with higher order symmetry 22 would have even greater 
load bearing capability. 

Discussion 

This paper paves the ways towards the study of interesting and 
unique geometric and mechanical properties of origami structures 
as mechanical metamaterials. It is expected that through a combina- 
tion of this approach and multiphysics simulations (e.g., COMOSL 
Multiphysics), more interesting properties can be explored. For 
example, the negative response between shear force and deforma- 
tion, and infinite tensile and bulk modulus may lead to some unique 
sound and vibration behaviors. When integrated with functional 
materials on origami patterns with micrometer feature sizes (e.g., 
the size of its rigid faces), such as nanowires and two-dimensional 
materials, the foldability of the origami pattern would provide 
unique tunable metamateirals with intriguing optical, electrical and 
magnetic properties, which is in fact under pursue. When combined 
with applications, the analysis of origami as mechanical metamater- 
ials can help to guide the development of origami based devices 8,9 . As 
all of these properties and applications are rooted from the way of 
folding, origami also provides a unique and powerful way on man- 
ufacturing. For example, plywood with unusually strong load bearing 
capability at the completed folded state can be manufactured in 
large-scale and low-cost by pre-creasing the wood panel based on 
the Ron Resch pattern. It is thus believed that origami may provide 
many interesting applications in science and engineering. 



Methods 

Non-local finite element method. Starting from the energy perspective, the elastic 
energy stored in a folded state is just the rotational energy at the creases since all the 
faces are rigid. If the creases are considered as elastic hinges, the elastic energy takes 
the quadratic term of the dihedral angels between creases. For the Miura-ori, the 
elastic energy can be written as 
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(e) 

Figure 3 | Load bearing of Ron Resch patterns, (a) A Ron Resch tube subjected to an axial compressive load, where the top view is given for the cross- 
section before the load is applied, (b) Normalized axial compressive force as a function of axial strain. Three representative states are shown as the insets at 
different strain levels. Their cross-sections and zoom-ins are also shown. Same scales are used in (a) and (b). (c) Illustration of a Ron Resch dome deforms 
to a completely collapsed state upon compressive load from the top, where the three-fold supporting structure is shown in the inset, (d) Photographic 
image showing the load bearing capability of a Ron Resch pattern at its completely collapsed state, (e) Photographic images showing the three-fold 
structures before (left panel) and after (right panel) the failure point is reached. The inset shows the instability, (f) Finite element simulation showing the 
first buckling mode of a Ron Resch plate with a three-fold supporting structure, (g) Finite element simulation showing the first buckling mode of a six- 
fold supporting structure. 
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where fc^ 0 and 0 are the spring constants of the hinges for dihedral angles oti and a 2 
for Miura-ori (superscript "Mo"), respectively; a 1>eg and a 2 , e9 are the corresponding 
dihedral angles for oti and a 2 at the undeformed state {or equivalently, just folded 
state); the summation runs over all dihedral angles. Similarly, the elastic energy can be 
readily constructed for the Ron Resch pattern, 



<=E5*f( 



EkH 



(8) 



where the superscript "RR" denotes the Ron Resch pattern and the subscripts have a 
similar meaning as explained for the Miura-ori. It is reasonable to take 
fej*» = ft*d = ft** and fef* = ^ = ^f = for paper folding (although in most 
machine-made papers, the fibers tend to run in one direction and so the hinge 
constants for edges running in different directions will be different). 

Because the dihedral angles are completely determined by the coordinates of 
vertices in rigid origami, the elastic energy can also be expressed as a function of 
coordinates of vertices, i.e., U tota i — U tota i(x), where x — (xi,X2, ■ ■ ■ ,x#) andx, is the 
position of a vertex i, and N is the total number of the vertices. When the external load 
F, is applied at vertex i, the total potential energy is Tl tgta i(x) — U tota i — Ffx,-. The 

equilibrium state of a rigid origami corresponds to a state of minimum energy and can 
be given by 



= o, 



(9) 



which needs to be solved to reach the equilibrium state of a rigid origami. There are 
many approaches that can be utilized to solve Supplementary Eq. S3, such as the 
conjugate gradient method and steepest descent method that just use the first deri- 
vatives of Yl totai , or the finite element method that uses both the first derivative (as the 

non- equilibrium force P= — — ^— ) and the second derivatives (as the stiffness 



matrix K = 



P l n lc 



dx 

- ) of n fof „;. The governing equation for the finite element method 



Ku = P, 



(10) 



where u = x — x (0) is the displacement of the vertices with x' 0) as the initial position of 
the vertices. For nonlinear systems, equation (10) is solved iteratively until the 
equilibrium characterized by the vanishing non- equilibrium force P — 0. For discrete 
vertices in rigid origami that has a great deal of similarity with atomic systems, the 
finite element method has been extended to capture the non-local interactions 23,24 . 

There are two aspects to consider when the finite element method is used. Firstly, to 
calculate the non- equilibrium force P and stiffness matrix K, the elastic energy U to tai 
needs to be explicitly written as a function of vertex coordinates, which is detailed in 
the Supplementary Information, Section "Nonlinearity of the Elastic Energy with 
respect to the Coordinates of Vertices". Therefore, iteration is needed. Secondly, non- 
local elements are required to capture the non-local interactions within a single 
element. For example, those nine vertices marked by blue, red and green dots in 
Fig. Id form one non-local element for the Miura-ori, focusing on the central vertex 
marked by the blue dot. Similar non-local elements {i.e., all solid circles and open 
circles in Supplementary Fig. S3) are constructed for the Ron Resch pattern, as shown 
in Supplementary Fig. S3. It may be noticed that the definition of non-local elements 
depends on specific rigid origami patterns and in each pattern different elements are 
formed for different types of vertices. These non-local elements are implemented in 
the commercial finite element package, ABAQUS, via its user defined elements 
(UEL). 
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